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ABSTRACT 



^ ' In this talk I will show that usual spin glasses are a peculiar kind of Abelian gauge 

O , theory. I will shortly review the techniques used to study them. At the end I 

^_( ' will consider more general models (e.g. spin glasses based on non Abelian gauge 

QQ , group) and I will discuss the relevance of these models to real glasses. Finally I 

^S| ' will derive from first principles a generalised Vogel-Fulcher law for the divergence 

of the characteristic time near the glass transition. 



1. Introduction 



The aim of my talk is to show how some of the subjects on which Oscar Klein 
t::^- I worked (in particular I refer to gauge theories) have been useful in apparently far 

^ ■ away subjects, i.e. spin glasses and real glasses. Indeed some models of spin glasses 

are a gauge theories of a rather peculiar type. 

I will start (in the second section) by briefly recalling the definition of a spin glass 
and presenting an idealised model. Later (in the third section) I will show that this 
^ ■ model is in reality a gauge theory corresponding to the Z2 group. Next (in the forth 

P ■ section) I will review the basic ingredients of the replica method which are used to 

study these models. I will later (in the fifth section) consider Abelian and non Abelian 
generalisations and I will derive from first principles a generalised Vogel-Fulcher law 
^ ■ for the divergence of the characteristic time near the glass transition. In the last 

section I will argue that these non Abelian models may share many properties with 
real glasses. 2. Spin Glasses 

I introduce here a very simple model of spin glasseditH. I consider a material where 
there are three kinds of atoms: M, A and B. M is a magnetic atom (it has a non 
zero magnetic moment) while A and B are magnetically inert. 

I suppose that at low temperature the system crystallises in such a way that the 
M-atoms stay on a regular lattice and the A and B atoms stay on the links of the 
same lattice. The position of the magnetically inert atoms is supposed to be random. 
In other words we consider an M A^ -Bioo-x alloy; the case a; = 50 corresponds to an 
equal proportion of A and B. 

Let us assume that the magnetic interaction among the M-atoms is of the nearest 
neighbour type and that it is mediated by the non magnetic atoms. The interaction 



among two M-atoms is ferromagnetic if the link is occupied by an A atom, while it is 
antiferromagnetic if the link is occupied hj a. B atom. We also assume for simplicity 
that the strengths of the ferromagnetic and of the antiferromagnetic interactions are 
equal. 

Usually the magnetic interaction is relevant only at temperatures much lower than 
the melting temperature and it may be neglected during the formation of the alloy. 
At low temperature the atoms may only oscillate around their equilibrium positions. 
If the temperature is decreased fast enough the position of the atoms in not influenced 
by the magnetic interaction. We can describe this situation by saying that we are in 
presence of a quenched disorder. 

If we assume that the spin are Ising variables, the corresponding Hamiltonian, in 
presence of a magnetic field h, is 

Hu[(r] = -^aiUi^k(^k-h^ai . (1) 

i,k 

The variables a are defined on the sites of the lattice and they take the values ±1. 
The variables Ui^k are defined on the links of the lattice, i.e. when i and k are nearest 
neighbours; they also take the values ±1. 

The variables U are random independent variables. For each choice of the U we 
can define a statistical expectation value: 

/^/ N\ Eaexp{-pHu[a])g{a) 
Eaexp{-(3Hu[a]) 

We are interested in computing the statistical expectation values averaged over the 
probability distribution of the samples, in other words the quantity 



{g{a))^{gia))u^J dPiU){gia))u, (3) 

where we denote by an horizontal bar the average over the U and P{U) is the proba- 
bility distribution of the variables U. We have already remarked that the probability 
distribution of the variables U does not depend from their interaction with the a: for 
each physical realization of the system they are fixed (quenched) variables. 

In the infinite volume limit the expectation value of intensive quantities does not 
depend on the realisation of the couplings U (i.e. all the samples have essentially the 
same properties) and the sample to sample fluctuations vanish in this limit. 

A very interesting quantity to evaluate is the magnetic susceptibility. A naive 
computation would give the following formula 



X = /?(1 - {a.)l) = P{1 - m{i)l) = P{1 - qea) , (4) 

where m{i)u = {(7i)u is the site dependent spontaneous magnetisation and qEA = 
m{i)'ij is the so called Edward Anderson order parameteilll. We shall see later how 
this formula for the susceptibility is modified by a more sophisticated treatment. 



Physical intuition tell us that at high temperature at zero magnetic field there is 
no spontaneous magnetization and consequently qea = 0. At low temperature each 
sample should develop its own spontaneous magnetization, and consequently qea 7^ 
at low temperature. The non vanishing of qea should therefore mark the spin glass 
transition. 

Before studying this model further it is convenient to analise its symmetries and 
in particular the consequences of gauge invariance. 

3. Gauge invariance 

Let us consider two different systems such that the couplings of the two systems 
are equal in all links not connected to a given point i and they differ by a sign only 
for the links connected to that point. These two sets of couplings have the same 
probability to be realised (all sets of couplings have the same probability). They have 
also the same free energy because the Hamiltonian does not change (at zero magnetic 
field) if simultaneously we change sign to the spin variable sitting at the point. 

In other words both the measures (on U and a) and the Hamiltonian are invariant 
with respect to the local gauge transformationa: 

o"i -^ -CTj ; Ui^k -^ Ui^k ■ (5) 

The set of all possible realizations of the system at zero magnetic field is gauge 
invariant under the gauge group Z2. The couplings and the spins play respectively 
the role of the gauge connection and of the matter field. At non-zero magnetic field 
the gauge invariance is explicitly broken. The Hamiltonian at zero magnetic field is 
(apart from a constant) the square of the covariant lattice-gradient of the a variables. 
It can be written as 

^^(CTi-f/,,,+^a,+^)2 . (6) 

i M 

The relevant quantities are gauge invariant. If two realizations of the couplings 
U and U' differ by a gauge transformation, their thermodynamic properties are the 
same. It is important to concentrate our attention on gauge invariant quantities. 

Let us study in more details how the thermodynamical quantities depend on the 
choice of the couplings. A quantity which is often used to characterise the gauge 
fields is the Wilson loop. We can associate to each closed circuit on the lattice the 
ordered product of all the links of the circuit. We thus define: 

W{C)= n U,,k. (7) 

{i,k)eC 

This quantity may take the values ±1. If W{C) = —1, the loop is said to be 
frustrates: along that loop it is not possible to find a configuration of spins such 
that 

Ui,kaiak = l y{t,k)eC. (8) 



If there are frustrated loops (i.e. if there is no gauge transformation which brings 
all couplings t/ to 1) it is not possible to find a configuration of the spins such that 
all terms in the Hamiltonian are positive. Some defects (i.e. links for which the 
contribution to the energy is negative) must be present. 

At low temperature the equilibrium probability distribution is concentrated on 
those spin configurations which have the minimal energy. It is interesting to study 
also those configurations which are local minima of the Hamiltonian in the sense that 
the Hamiltonian increases when we flip a spin. These local minima are very important 
in the dynamics outside equilibrium because at low temperatures the system may be 
trapped for a very long time in these minima. 

If we study the structure of local and global minima with great care, we discover 
that frustration implies the presence of defects which can be put in many ways on the 
lattice. The ground state is degenerate. The number of local and not global minima 
is also large. 

This phenomenon is well known in gauge field. On the lattice the choice of the 
Landau gauge corresponds to find the maximum of 

Y^Tig^U.^ugk. (9) 

i.k 

where g is the gauge transform which brings the gauge fields in the Landau gaugea. 
This problem is equivalent (in the our case) to find the minimum of the Hamiltonian 
eq. (|1|). Gribov ambiguity tell us that in the general case there are many possibility 
of choosing the gauge. This result implies the existence of many minima of the 
Hamiltonian eq. (|I|)f 

4. The Replica Method 

We face now the problem of evaluating the quantities which appear in equation 
(0). The first proposal would be to sum over the variables U and to remain with an 
effective interaction for the variables a. This approach clashes with fact that both 
the numerator and the denominator of eq. (|^) depend on the variables U and the 
sum over the U is not easy. 

This difiicultv may be avoided by introducing n identical copies (or replicas) of 
the same system El. We define 

,(^\ Eg Ec7 exp(-/?//n[ff, U])g{a^) , . 

L<7 L(7 exp(-/yii„[o-, u\) 

where the spins af carry an other index (a) which ranges from one to n. The new 
Hamiltonian is the sum of n identical Hamiltonians 

Hr,ia,U)= J2 HuW"]- (11) 

a=l,n 



*In the continuum Gribov ambiguity is normally present in non Abelian gauge theories. It is also 
present in Abelian theories when we allow configurations which are singular in the continuum limit, 
e.g. like magnetic monopoles. 



It easy to check that 



w.)).. ^ ^"<f'>if'" . (12) 



where the [/-dependent partition function is defined as 

Zu = Y.^xp{-(3HuW]). (13) 

O" 

We finally find that 

{g{a)) = {g{a))n\n=o ■ (14) 

In this way one finds that properties of the matter fields averaged over the dis- 
ordered (quenched) gauge fields can be computed by considering the gauge fields 
interacting with n copies of the matter field and computing the expectation values in 
the limit n — *> 0. The argument is familiar to those who work in numerical simulation 
of latttice gauge theories, where n plays the role of the number of quark fiavours in 
the sea! 

The average of the U fields can be done and one remains with an effective inter- 
action for the a variables. This effective interaction must be written in terms of the 
gauge invariant combinations. In this case the most appropriate variables are 

Qf = <^l (15) 

The diagonal terms of the matrix Q are identically equal to 1, so that we can consider 
only the off-diagonal term. All computations must be done for generic values of n 
and we must send n to zero at the end. 

As usually we can expand the effective interaction in powers of Q, being careful 
to preserve the various symmetries of the problem, i.e.: 

• The group of permutations of the n replicas S{n). 

• The spin reversal symmetry for each replica, i.e. n times the direct product of 
the Z2 global group, where each Z2 group acts on a different replical 

In the continuum limit in D dimensions the simplest form for the effective free 
energy, in which all important terms are contained, is the following : 

F[Q] = jd^x [^Y.^T{d,Q{x)f + W{Q{x))\ , (16) 



M 



where the function W{Q)is given by 



W{Q) = rTr(Q2) + g,Tr{Q') + g,Ti{Q') + yY.Qli> ■ (17) 

ab 



^It is evident why I proposed the name quenched for the approximation of neglecting quark loops in 

QCD. 

■I- This symmetry is present only at zero magnetic field. 



The usual strategy consists in computing the minimum of F and constructing the 
usual perturbative expansion for the small fluctuations around this minimum (the so 
called loop expansion). As we shall see even the first step is not very simple. 

Let us assume that the F[Q\ is minimized by a function Q{x) which does not 
depend from x. We can this set Q{x) = Q and look for the minimum of W{Q). 

For r > we find that there is only a minimum at Q = 0. When r is negative we 
easily find that there is a stationary points at 

Qab = q, y ab. (18) 

This point is invariant under the replica group S{n). In order to decide if this sta- 
tionary point is a minimum of W, we have to compute the small fluctuations around 
this point. The Hessian 

dW 

'Hah,cd = r,^ ^ , (19) 

must have non negative eigenvalues (in a field theory language there should be no 
negative squared masses). 

If we consider the case where y = 0, the S{n) symmetry is promoted to an 
0{n) symmetry (i.e. there is an accidental symmetry). This 0{n) symmetry is 
spontaneously broken by the choice in eq. (|18|). Consequently the Hessian has zero 
eigenvalues (i.e. there are Goldstone Bosons). 

If y is different from zero, the 0{n) symmetry is explicitly broken and the Gold- 
stone Bosons acquire a mass squared proportional to y. A detail computation shows 
that in the general case y is negative and the Goldstone Bosons acquire a negative 
mass squared (i.e. the Hessian has negative eigenvalues). 

This instability implies that the proposed stationary point is not a local minimum 
and one has to look for a minimum which williie no more invariant under the Sn 
group. Such a minimum can be constructed H'Q; we have to introduce an infinite 
sequence of steps breakings the S{n) symmetry. After a rather long computation one 
finally finds a matrix Q to which one can associate a function q{u), where u belongs 
to the interval — 1. 

There is a simple physical interpretation of this symmetry breaking. The presence 
of frustration implies that for each realization of the couplings there are different equi- 
librium states with different local magnetisation (i.e. different vacua). We indicate 
by w" the probability of finding the system in the state labelled by a and by mf 
the local magnetization (i.e. the expectation value of ctj in the state a). The overlap 
among different states may be defined as 

g-^ = ^^=i.^<^^ (20) 

After some computations one finds that 

^«;"^7/(g"7) ^ j dqP{q)f{q) = [' dufiqiu)). (21) 



The probability of having two states with given overlap q is thus controlled by the 
function q{u). 

This approach is able to explain in a qualitative and sometimes a quantitative 
way many of the properties of real spin glass. When the replica symmetry is broken, 
a change in the external magnetic field changes the structure of equilibrium states 
and one must be quite carefully in the definition of the susceptibility. One finds that 
there are two susceptibilities: 

• The linear response susceptibility which quantifies the response to a small varia- 
tion of the magnetic field measured on a short time scale such that no global re- 
arrangement of spins is possible. This susceptibility is given by xlr = P{^ — <1ea- 

• The thermodynamic susceptibility which quantifies the response to a small vari- 
ation of the magnetic field on a very long time scale such that global rearrange- 
ments of spins are possible. This susceptibility is larger^ and it is given by 
Xeq = /3{1- fduq{x)). 

The difference between these two susceptibilities is one of the most characteristic 
phenomena experimentally observed in spin glasses and it is well explained by this 
theoryB. 

At this stage of the theory, where fluctuations are neglected, is not clear which of 
the theoretical predictions does survive when the effect of the fluctuations is included. 
In order to have more precise theoretical prediction for real three dimensional systems 
it is necessary to go beyond the mean field approximation. This has been the subject 
of intensive studies. Unfortunately the computation of the corrections to the mean 
field theory are rather involved. The correlation functions at zero loops are rather 
complicated Q and the one loop corrections to the correlation functions have not yet 
fully computed. 

Numerical simulations in 4 dimensions Bi strongly suggest the correctness of the 
broken replica picture. No anomaly is observed. Numerical simulations in 3 have not 
produced a compleately clear picture, although some progresses have recently been 
done 0§. 

When we will master the loops corrections, we will be hopefully able to set up a 
renormalization group study of the properties of the system in the low temperature 
region and in particular to compute the value of the lowest critical dimension, i.e. 
the dimension where the structure of replica symmetry breaking scheme is no more 
consistent and the predictions of the mean field theory do not apply anymore?' 

5. Other Gauge Groups 

There is no need to consider only the Z2 group. The generalisations to other 
groups are very interesting. In particular the case of the U(l) group is relevant for 

§Onc finds that qea ~ max(g(u))). 

^For the 0{n) symmetry the lowest critical dimension is 2, for the Z2 symmetry it is 1. 



the study of irreversibility in high Tc superconductors 113. Here we consider a different 
gauge group, non-Abehan, for which new results can be found. 

In this model each spin can takeiivaluesi^l. The gauge field Ui^k is an element of 
the permutation group of p elementsta. The Hamiltonian is 

Hu[<^]=T.^<r.u.,.,- (22) 

i,k 

The gauge group is the group of permutation of p objects (S'(p)). If p = 2 we recover 
the previous case (Z2 = >S'(2)). 

Due to the non Abelian nature of the gauge group there is no global S{p) symmetry 
at fixed U . As an effect of this decrease in the symmetry the effective free energy 
contains an extra cubic term 

W{Q) = TTr{Q') + gsTr{Q') + zJ2qI, + 0{Q'), (23) 

ab 

which is forbidden in the Abelian case at zero magnetic field. 

The value of z is crucial. If z is small with respect to g^ {z = at p = 2) we find 
that the value of Q at the saddle point is a continuous function of the temperature. 
For large values of z (in mean field we need 2; > 4) we fins that the value of Q jumps 
discontinuously at the phase transition T^. 

A dynamical computation shows that in these conditions there are metastable 
states above the critical temperature, which correspond to vacua having high free 
energy of the true ground states. The system rnay be trapped in these false vacua, 
and it may remain there for a quite long time Elll3. 

If we consider only point independent minima of the effective Hamiltonian we find 
that there is a dynamical temperature (To) at which the characteristic time (tc) of 
the system divergeslll. 

Below Tu the characteristic time diverges exponentially with the size of the system. 
This exponentially large time may be computed as follows. We introduce an effective 
potential V{q) defined as: 

- NViq) = I dP[a] In (| dP[fi]6iq - g.,j) , (24) 

where 

i^^f^ = j;^ (25) 

and for simplicity we use the short notation 

[dP[a]=Y.exp{-(3Hu[a]). (26) 



Using the techniques of llaO this potential can be computed by considering the 
effect of coupling R = 1 + e replicas in the limit e — > 0. It turns out that this potential 
has always a minimum at g = 0. For T < To we finds that this potential acquires a 

8 



new minimum at g = g^, the two minima being obviously separated by a maximum 
at g = Qm- 

If we assume that the dynamics may be approximated by the motion of the system 
in this potential, we find that 

t,^exp{N6{T)), (27) 

where the function 6{T) is different from zero below To and is given by 

SiT) = ViqM)-V{qm), (28) 

i.e. by the free energy barrier that the system has to cross. 

A characteristic time which increases exponentially with the size of the system 
cannot be the correct answer in a system with a short range Hamiltonian. It is 
well known that in this case metastable states do not exist: false vacua do decay by 
thermodynamic tunnelling and the characteristic time can be computed by evaluating 
the free energy of the critical droplet (or instantons in a field theoretical language). 

If we use the thecnique of c2l to study the instantons of these theories we find out 
that the characteristic time is finite also below Tq: 

t,~exp(A(r)), (29) 

where the function A(T) diverges at the critical temperature as 

A(r) oc (T - Te)-(^"i) (30) 

These results can be qualitatively easily explained. The situation is quite similar 
to the usual enucleation theory in which we have to estimate the size of the critical 
droplets. The time for escaping from the false vacuum is given by the exponential of 
the instanton action (or the exponential of beta times the free energy barrier). An 
instanton of radius R has two contributions to its action 

A{R) = -tR'^ + Si?(^-i) (31) 

The first is a bulk contribution, which vanishes at the critical temperature (r oc T—T^) 
and the second is an surface term (the free energy of the interface). The radius for 
which A{R) is maximum diverges as the inverse of r and the free energy excess at 
this point is proportional to (T — Tc)~^^~^\ 

Near the critical temperature the evolution of the systems is dominated by the 
rearrangement of spins in large domain, whose radius R diverges at the critical tem- 
perature. Larger and larger barriers must be crossed when we go near T^ and the 
characteristic times diverges at this temperature. 

In the Z2 case discussed in the previous section there are no metastable states 
already in the mean field approximation c3t3. In short range models the characteristic 
time is divergent at the transition, but it divergent as a power law, i.e. as {T — Tc)~''^, 
where A = 2 in the mean field approximation. The new feature of the non Abelian 

9 



models is the exponential divergence of the characteristic time at the phase transition. 
This divergence makes extremely hard to thermalize the systems also at temperatures 
not so close to T^. 

6. Toward Real Glasses 

In spin glasses the Hamiltonian is random as an effect of quenched random dis- 
order. In real glasses the Hamiltonian is not random and the quenched disorder 
is dynamically generated at low temperature. We can therefore ask how much of 
the qualitative and quantilatite results which have obtained in spin glasses may be 
transferred to real glasses E3~^. 

In order to understand this point one can study models in which the Hamiltonian 
does not contain quenched disorder and to compare the results with those coming of 
random HamiltonianE3~Ea. 

In general our strategy is the following. We want to study the properties of a 
given Hamiltonian Hg which is not random. We consider a class of Hamiltonians 
Hfi, of which Hg is a particular case. We choose the class H^ in such a way that the 
statistical properties of Hg and that of a generic Hamiltonian in Hpi are as similar 
as possible. In the best case we can obtain that the two corresponding free energies 
coincide in the high temperature expansion. In general the behaviour of the system 
can be controlled better in the high temperature phase. 

After having constructed H^ in an appropriate way, we can suppose that the 
statistical properties of Hg and Hji are the same or, if they are different, we can 
construct a perturbative expansion which compute this difference. It is clear that 
this approach may be successful in the high temperature region (more or less by 
construction) and it may also reproduce the behaviour of the system in the glassy 
region, included the dynamic and static transitions. However it is cannot certainly 
reproduce the possible existence of an ordered crystal phase. 

The proposed strategy for approximating a given Hamiltonian with a random 
Hamiltonian works very ^^llin the many cases, at least in the framework of the 



mean field approximationE3~E3. Real glasses have a qualitative behaviour that is 
quite similar to the one described in the previous sectionP It is therefore tempting to 
suppose that also real glasses belongs to the same universality classes of system with 
a random Hamiltonian. If this conjecture is correct equation (^) should apply also 
to real glasses. 

It well known that glasses show an extremely large increase in the value of the 
viscosity (which is proportional to the characteristic time) near the glass transition. 
This increase of fifteen order of magnitude has been measured in many materials. 
In the case of fragile glasses this increase has been fitted by the phenomenological 
Vogel-Fulcher law E3: 

r/ocexp(C(T-r,)'i). (32) 

II Theoretical arguments can be done which point in this direction l3. 
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The prediction of the rephca approach for a three dimensional systems is 

77cxexp(C(r-re)-2). (33) 

Unfortunately the data are not sensitive enough to make a clear cut choice among the 
two proposed formulae, although the data are fitted slightly better by eq. (|33D than 
by the traditional Vogel-Fulcher law. More careful studies of this point are needed ** 

It is quite amazing that starting from the study of some peculiar forms of gauge 
theory on the lattice we finally end up with a partial theoretical understanding of the 
phenomenological Vogel-Fulcher law, which have been already proposed when Oscar 
Klein was young. 

The possibility of applying these geometrical ideas (which originate from particle 
physics) to such far away systems as spin glasses and real glasses, shows the substantial 
unity of theoretical physics, a field to which Oscar Klein gave so many and important 
contributions. 
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